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ABSTRACT 

A comparison between published field galaxy stellar mass functions (GSMFs) shows that 
the cosmic stellar mass density is in the range 4-8 per cent of the baryon density (assuming 
£!b = 0.045). There remain significant sources of uncertainty for the dust correction and un- 
derlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass 
function. We determine the z < 0.05 GSMF using the New York University - Value-Added 
Galaxy Catalog sample of 49968 galaxies derived from the Sloan Digital Sky Survey and var- 
ious estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is 
fitted with a double Schechter function that has a>2 ~ —1.6. At masses below ~ 1O 8 5 A40, 
the GSMF may be significantly incomplete because of missing low surface-brightness galax- 
ies. One interpretation of the stellar mass-metallicity relation is that it is primarily caused by 
a lower fraction of available baryons converted to stars in low-mass galaxies. Using this prin- 
cipal, we determine a simple relationship between baryonic mass and stellar mass and present 
an 'implied baryonic mass function'. This function has a faint-end slope, 012 — —1-9. Thus, 
we find evidence that the slope of the low-mass end of the galaxy mass function could plausi- 
bly be as steep as the halo mass function. We illustrate the relationship between halo baryonic 
mass function — » galaxy baryonic mass function — > GSMF. This demonstrates the require- 
ment for peak galaxy formation efficiency at baryonic masses ~ lO n .A40 corresponding to 
a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower 
masses. 

Key words: galaxies: evolution — galaxies: fundamental parameters — galaxies: halos — 
galaxies: luminosity function, mass function. 



1 INTRODUCTION 

The galaxy luminosity or mass function is a fundamental tool 
used in interpreting the evolution of galaxies. The functions are 
usually defined as the number density of galaxies per logarith- 
mic luminosity or mass interval. A steeply rising mass function 
to the faint population has been a generic prediction of galaxy 
formation based on hierarchical clusterin g iWhite & ReeslTl97a : 
iKauffmann elaH 1 19931 ; ICole et alj 1 1994) . In contrast, the field 
galaxy luminosit y function was observed to have a signi ficantly 
flatter 'faint end' dBinggeli et al.ll988l ; lLovedav et aTTl 992). In or- 
der to reconcile cold dark matter (CDM) galaxy formation mod- 
els with the observed luminosity function, star formation (SF) 
is suppre ssed in low-mass halos by, for ex ample, supernovae 
feedback iLacev & SilkHl99ll;lKav et alj|2002h or photoionisation 
fcfstafhioiil 1992l ; ISomervilld2002l) . ~ 

On the observational side, accurately determining the num- 
ber densities of faint galaxies in the field or in clusters is challeng- 
ing. This is mainly because of the low surface brightness (SB) of 
these galaxies, which means that they can be undetected in pho- 



tometry even if h igh-SB galax i es with the same apparen t magni- 
tude are detected ( lDisnevlll976] ; lDisnev & Phillipp Jl983h . Despite 
these and other challenges, there is evidence for a 'faint-end up- 
turn' in luminosity functions whereby the luminosity function is 
rising steeply fainter than about 3-5 magnitudes below the char- 



acteri stic luminosi ty, in clusters Driver et aU 



2005|) or the field jLovedaylll997l ; IZucca et al 



19981 : iBlanton et aflT2005a[). Note however that the upturn is n ot 
always evident (Norberg et al. 2002; Harsono & de Pro pris 2007). 



19941: JPopesso et al 



1997l : lMarzkeetal 



A faint-end upturn suggests that the efficiency of feedback 
has levelled off at low galaxy luminosities. A more direct analy- 
sis is to compare galaxy mass functions with predicted mass func- 
tions from CDM models. iKlypin et alj J1999 ) compared the circu- 
lar velocity distribution of satellite galaxies in the Local Group (see 
also lMoore et alj 19991) . This highlights the 'substructure problem' 
where there are 5-10 times as many low-mass satellites in CDM 
models than observed. However, we still have not reached a com- 
plete census of Local Group satellites as evidenced by recent dis- 
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coveries of low -SB galaxies around M 31 dMartin et al]|2006h and 
the Milky Way dBelokurov et al.ll2007b . 

Reliable dynamical mass estimates for complete and large 
samples of field galaxies are difficult to obtain. Stellar masses 
can be estimated for significantly larger sample s of galaxies us- 
ing the principles o f stellar population synthesis dTinslev & Gunnl 
1 1 976L iTinslevfl 1 98C>) . Thus, galaxy stellar mass functions (GSMFs) 
can be derived from galaxy luminosities (e .g., Balog h~et alj|200ll ; 
ICole et al.ll200ll ; lBell et afll2003bl : lKodama & Bowej|2003l) . Com- 
parisons between GSMFs are then, in theory, free of the stellar- 
population contribution to mass-to-light (M-L) variations that are 
inherent in comparisons between luminosity functions. 

Integrating a field GSMF to determine the cosmic stellar mass 
density (SMD) brings to light the large gap between this value 
(fistars ~ 0.003) a nd estimates of Q.b from Big-Ba ng Nucleosyn- 
thesis theory (0.04, iBurles, Nollett. & Turned [200 1 ) or the power 
spect rum of the cosmic background radiation (0,045, [Spergel et all 
120071) . Even determination of the baryonic content of galaxies in- 
cluding stare_aj^ coldgas accounts for less than or about 0.1 
of fi 6 jBelletalJ 120033) . About 0.2 of Qb is accounted for by 
hot plasma identified by X-ray emission in clusters and groups 
(Fukugita, Hogan, & Peebles 1998); while the rest is projected to 
be in a more diffuse inter-galactic medium (Cen & Ostriker 1999) 
and perhaps the 'coronae' of galaxies ((Mailer & Bu llockll2004l) . 
The overall efficiency of baryons falling into the luminous disks 
or bulges of galaxies is low. 

The efficiency of SF (fraction of baryonic mass converted to 
stars) varies significantly with galaxy mass. A low efficiency aver- 
aged over the life of a galaxy gives rise to a low gas-phase metal- 
licity (in the inter-stellar medium) because the metal pro duction 
by supernovae is diluted by gas reservoirs within a galaxy (Tinslev 
ll980luBrooks et al.l2007l) and further infalling material. Conversely 
a high efficiency can drive the gas-phase metallicity to high values. 
In this paper, we explore th e use of the stellar mass-metallicity re- 
lation ( iTremonti et alj2 004) as a SF efficiency estimator to convert 
the GSMF to a baryonic mass function. 

The plan of this paper follows. Determinations of the field 
GSMF and the SMD are reviewed in §|2] An estimate of the low- 
redshift field GSMF paying careful attention to SB selection effects 
is described in §[3] The evident faint-end upturn is compared with 
cluster GSMFs. The relationship between gas-phase metallicity and 
stellar mass is then used to convert the field GSMF to an 'implied 
baryonic mass function' assuming a simple relationship between 
metallicity and the fraction of baryonic mass in stars. This is de- 
scribed in §|4] This galaxy baryonic mass function, and the GSMF, 
are compared with halo mass functions and related to galaxy for- 
mation efficiency as a function of mass. Summary and conclusions 
are presented in § [5] The dependencies of stellar M-L ratio on var- 
ious assumptions are presented in the Appendix. Throughout this 
paper, a cosmology with Ho = 70kms _1 Mpc -1 , fl m ,o = 0.3 
and Qa.o ~ 0.7 is assumed. 



2 COMPARISON BETWEEN PUBLISHED GALAXY 
STELLAR MASS FUNCTIONS 

By estimating stellar M-L ratios for galaxies it is possible to calcu- 
late the equivalent of a luminosity function, for the total s tellar mass 
of ga laxies, known as the galaxy stellar mass function dCole et all 
1200 lh . This then gives a more fundamental account of the baryons 
that are locked up in stars and how they are distributed amongst 
galaxies of different masses. 
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■ Colo et al. 2001 table 4 
Bell et al. 2003 table 4 

■ Bell et al. 2003 table 5 
Baldry et al. 2004 fig. 10 

■ Baldry et al. 2006 fig. 8 
Panter et al. 2007 fig. 3 
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Figure 1. Comparison between published field galaxy stellar mass func- 
tions (z < 0.1). The mass functions have been conv erted to Ho = 
70kms _1 Mpc -1 where necessary and the ICole et alj GSMF has been 
converted to the diet Salpeter IMF. 



One of the key ingredients in this calculation is the stellar ini- 
tial mass function (IMF), which is generally assu med to be inde- 
pendent of galaxy type or mass ( El me green 1 200 lh . For the com- 
parisons in this paper, the IMFs used are the 'diet Salpeter' that 
is defined as . 7 x the mass derived from a stand a rd Sa lpeter 
l lBell & de Jon3l200ll : iBell et alJl2003bT) . the lKroupal feOOll) . and 
the lChabriea ( 120031) IMF. These are all similar in terms of M-L ra- 
tio as a function of galaxy colour. The variations highlighted and 
discussed between different mass estimates in this paper are not 
significantly dependent on IMF choice (Appendix). 

Figure[T]shows a comparison between published field GSMFs 
(where 'field' in this case means a cosmic volume average). A 
brief description of how these were derived is given below. Masses 
and number densities were converted to a cosmology with Ho = 
70 km s _1 Mpc -1 where necessary. 

(i) IColeetalJjjxiil) : data were taken from the Salpeter column 
of table 4. The masses were multiplied by 0.7 to convert to the diet 
Salpeter IMF. The published GSMF was derived from a match of 
the T wo Micron All Sk y Survey (2MASS) extended source cata- 
logue Jj^t^^lHoOa) to the 2dF G alaxy Redshift Survey (2dF- 
GRS: IColless et al.l200ll) . ICole et al.1 computed bj - K and J - K 
colours for a range of exponenti ally declining SF h i storie s with 
Zform = 20 and metallicities using lBruzual & Chariot ( 1993) mod- 
els. Dust attenuation from the models of iFerrara et alj d 19991) were 
applied. Stellar M-L ratios in the near-IR bands were determined 



from the models that most closely matched the observed galaxies' 
colours; with stellar masses assumed to be 0.72 times the integral 
of the SF rate (to accoun t for recycling). 

(ii) IBell et al.1 d2003bl) : data were taken from table 4. The plot- 
ted GSMF is the sum of the g-band derived early- and late-type 
Schechter functions. The GSMFs were derived from a match of 
the 2MASS catalogue to the Sloan D igital Sky Survey (SDSS; 
I York et al.1 [20001 : IStoughton et alj|2002h . The ugrizK data were 
fitted to magnitudes computed for a ra nge of exponential SF his- 
tories and m etallicities using PEGASE dFioc & Rocca-Volmerangel 
1 19971 ll 999) models. M-L ratios were determined from the best-fit 
models. 

(iii) IBell et ail as above but using data taken from table 5. This 
is a bin ned GSMF d e rived from the same M-L ratios. 

(iv) iBaldrv et al.1 d2004l) : data were taken from fig. 10. The 
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GSMF is the sum of the red- and blue-sequence Schechter func- 
tions. This was derived from SDSS data with a M-L ratio given by 
a linear function of colour: log(A4/Lr) = —0.55 + 0.45 (it— r) U 
This relat i on was det ermined from stellar m asses computed by 
iBell et al.l J2003bh andlKauffmann et ail d2003h . 

(v) iBaldrv etd] J2006I) : data were taken from fig. 8. The GSMF 
is the sum over all the environments . This was derived in a simi- 
lar way to that of Baldrv et al. ( 2004) using a different relationship 
between M-L ratio and colour. In this cas e, the relation wa s deter - 
mined from stellar masses computed by Kauffmann et alj J2003t) 
and K. G. These me t hods a re described in §[3] 

(vi) IPanter et al.l fc007l) : data were taken from fig. 3. The GSMF 
is derived from the Schechter parameters. The method for estimat- 
ing stellar masses is described in §[3] 

The GSMFs are similar except for the IBaldrv et all fc006h 
function, which appears to be shifted by about 0.2 dex to lower 
masses. This is because the method used by K.G. for that paper un- 
derestimated the masses of luminous red galaxies. This may reflect 
degeneracies with fitting to only ugriz photometry and attempting 
to fit SF bursts, metallicity and dust attenuation (Appendix). 



2.1 The stellar mass density of the Universe 

A quantitative overall comparison can be made by integrating each 
GSMF to obtain the total SMD. The critical density for H = 
TOkms" 1 Mpc" 1 is p c = 1.36 x K P.MpMpc -3 and we as- 
sume n b = 0.045 JSpergel et alj|2007h . The fraction of baryons in 
stars f s = fi s tarsM is then (i) 0.054, (ii) 0.061, (iii) 0.063, (iv) 
0.056, (v), 0.035, (vi) 0.042. Note these values do not include ex- 
trapolations of the GSMF to lower masses than the plotted points or 
corrections to total galaxy magnitudes from catalogue magnitudes 
unless already applied^ 

The total SMD can also be determined usi ng methods 
other than integrating a GSMF. iDriver et all ( 120071) determined 
f a = 0.083 ± 0.012 from the B-band luminosity function us- 
ing an empirical geometrica l-based correction for dust attenuation; 
Baldr v & Glazebrookl J2003h give a range from 0.035 to 0.063 from 
fitting to the local cosmic luminosity densities (over a range of 
IMFs but similar to th ose used here in terms of M-L ratio); and 
iNagamine et al.U2006l) found 0.051 ± 0.009 from fitting to various 
d ata including the extr agalactic background light. See also table 5 
of lGallazzi et alj d2008l) for a list of SMD measurements. 

Overall there remains a considerable uncertainty in f2 s tars or 
f s , which is probably in the range 4-8 per cent (for the Chabrier, 
Kroupa or diet Salpeter IMF). This is much larger than the uncer- 
tainty in Q m or Qt within the standard cosmological paradigm. 

Accounting for dust is obviously a major problem. Correct- 
ing for attenuation by incorporating a dust law in stellar popula- 
tion synthesis (PS) models, and fitting to colours and/or spectral 
features, cannot account fo r populations behind opaque screens. 
Notably. IDriver et"al] 12007) obtain the largest value of f s and the 
method uses comparisons of disks and bulges viewed at v arious ori- 
entations to constrain their dust model ( Tuffs et all2004l) . It is then 
useful to compare the SMD with near-IR measurements that are 



1 M or M s is used for stellar mass in this paper. 

2 SDSS Petrosian magnitudes theoretically recover 99% of the flux of a 
galaxy with an exponential profile, 80% with a de-Vaucouleurs profile, and 
95% in the case of nearly un resolved systems (c l ose to the point-sprea d 
function of SPSS imaging) <Blanton et ai]|200ll ; Istoughton et alj|2002h . 
lBelletalj f2003b) have used estimated corrections for this loss of light. 



minimally affected by dust attenuation. The total A'-band luminos- 
ity de nsity, in units of 10 8 I/p)Mp c~ 3 , is given by 4 .0 I Cole et all 
l200lh . 5.0 jKochanek et al]|200ll) and 4.1 IBell et alj|2003bl) l 3 l Us- 
ing these values and the full range of f s = 0.035-0.083 implies a 
cosmic \og{M/LK) in the range —0.35 to 0.1. A high SMD with 
fs ~ 0.08 appears unlikely given that only the oldest simple stellar 
populations have \og(Ai / Lk) — 0.0 (Appendix); unless there is 
more than nominal attenuation of the A'-band and/or the A'-band 
luminosity density is underesti mated. 

How can we reconcile the lDriver et al. | i[2007j) result (0.083 ± 
0.012) with the standard PS methods of Fig.Q](0.04-0.06 ignoring 
extremes) other than resorting to saying the higher result is 2-sigma 
too high? The standard PS methods cannot account for populations 
be hind optically -thick screens and use less deep photometry than 
the IDriver et alj result, which is based on the Millennium Galaxy 
Catalogue (MGC). A reasonable estimate of unaccounted for light 
would be 20% bringing f s up to 0.05-0.07 for the standard PS 
methods0On the other hand, the MGC estimate is subject to larger 
cosmic variance and the luminosity density appears to be high by 
lO-2O%0and the M-L ratios assumed for the attenuation-corrected 
magnitudes could also be high by 10%. Accounting for these brings 
fs down to 0.06-0.07 in agreement with the standard PS methods. 

While the largest contribution to the SMD comes from galax- 
ies around the break in the GSMF, lower mass galaxies play a key 
role in the processing of baryons. In the next section, we discuss a 
new calculation of GSMFs using various stellar mass estimates and 
consider how accurately the lower-mass end can be determined. 



3 GALAXY STELLAR MASS FUNCTIONS FROM THE 
NYU-VAGC LOW REDSHIFT GALAXY SAMPLE 

A large low-redshift sample derived from the SDSS is the New 
York University - V alue-Added Galaxy Catalog (NYU-VAGC) 
( iBlanton et"ai]|2005allbl) Fl While the data are obtained from stan- 
dard SDSS catalogues, the images have been carefully checked 
for artifacts including deblending. The data cover cosmological 
redshifts from 0.0033 to 0.05 where the redshifts have been cor- 
rected for peculiar velocities using a local Hubble-flow model 
dWillicketal]|l997h . We use the NYU-VAGC low-z galaxy sam- 
ple to recompute the GSMF down to low masses. 

The Data Release Four (DR4) version of the NYU-VAGC 
low-z sample includes data for 49968 galaxies. These da t a wer e 
matched to stellar m asses estimated by iKauffmann et al.1 d2003l) , 
iGallazzi et al] fc005h . and lPanter etafl J2007h ; with 49473, 32473 
and 38526 matches, respectively. Minor adjustments were deter- 
mined to account for variations between the SDSS data from which 



3 iHuang et alj J2003h reported a high value for the ft'-band luminosity 
density of ~ 8.0 (lO 8 L0Mpc~ 3 ). However a reanalysis by K.G., of the 
Hawaii-AAO data extended to K < 16, optimised for the accuracy of the 
luminosity density at z < 0.2 gives a value of ~ 4.5 in agreement with the 
2MASS results. 

4 The estimated corrections to the SMD values derived from the GSMFs of 
Fig-Hare about 5-10% to account for missing low-SB galaxies, low-mass 
galaxies and/or corrections to total magnitudes, and 10-15% for optically- 
thick regions. The latt er corresponds to ap proximately the attenuation in the 
i4"-band derived from lDriver et "aI U2007t) . 

5 The uncorrected total B-band luminosi ty densities are 1.6 from MGC 
iDriver et all2007l) and 1.3 from 2dFGRS jNorberg et alll2002l) in units of 
10 8 L Q Mpc- 3 (H = 70kms _1 Mpc" 1 ). 

6 NYU-VAGC data are available from 
http://sdss.physics.nyu.edu/vagc/ 
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the stellar masses were determined and the NYU-VAGC (photom- 
etry and cosmological redshifts). Where no stellar mass was avail- 
able for a galaxy, the stellar mass was determined using a colour- 
M-L relation calibrated to the particular set of stellar masses. In 
addition, stellar masses were computed by K.G. using the NYU- 
VAGC Petrosian magnitudes. Thus there are four stellar masses for 
each galaxy. The methods are described briefly below. 

(i) Kauffmann: The stellar masses were obtained by fitting a 
grid of population synthesis models, including bursts, to the spec- 
tral features D4000 and H<5 absorption. The predicted colours were 
then compared with broad-band photometry to estimate dust atten- 
uation. Stellar M-L ratios were determi ned and applied to the Pe t- 
rosian z-band magnitude. For details see Kauffmann et al. (2003)0 

(ii) Gallazzi: The computation was similar to the above method 
except five spectral features were used. The features were carefully 
chosen for calculating stellar metallicity and light-weighted ages 
while mini mising effects caused by chemical abundance ratios. For 
details see lGallazzi et al . (2005)Ef 

(iii) Panter: The computation involves fitting synthetic stellar 
populations to eac h galaxy spectrum using the MOPE D data com- 
pression technique l lHeavens. Jimenez. & Lahavl200ri) . The SF his- 
tory of each galaxy was modelled using 1 1 logarithmically-spaced 
bins in time, each with an SF rate and metallicity, and a sim- 
ple dust screen. The tota l stellar mass is that obtained from the 
iBruzual & Chariot! (2003) models given the best- fitting SF history 
as input. For details see iPanter et al. N2004l2007l) and references 
therein. 

(iv) Glazebrook: This is the only purely photometric method. 
Stellar masses were determined by fitting to the observed-frame 
Petrosian ugriz magnitudes of each galaxy. A grid of colours were 
computed using PEGASE models. A range of exponential SF his- 
tories and metallicities were input. Bursts were added with mass 
ranging from 10~ 4 to t wice the mass of the primary component 
( Glaze brook e7ai]|2004l) . For this paper, no dust attenuation was 
included in the models (see Appendix for discussion). Depending 
on the attenuation law, incorporati ng dust can be close t o neutral in 
terms of M-L ratio versus colour (Bel l&deJongll200ll ; see fig. 12 
of iDriver et al 1 l2007l for the effect modelled from face-on to edge- 
on). 

The IMFs used were the similar Kroupa and Chabrier IMFs, 
for methods (i,iv) and methods (ii,iii), respectively. Comparing the 
different mass estimates for each galaxy, the standard deviation is 
typically in the range 0.05-0. 15 dex. The standard deviation is gen- 
erally lower for red-sequence galaxies compared to blue-sequence 
galaxies. 

For each set of stellar masses, the galaxies were divided into 
logarithmic stellar mass bins. For each bin, the GSMF is then given 
by 



(1) 



Alog.M ^ V m 



where: Vmax.i is the comoving volume over which the ith galaxy 
could be observed, and Wi is any weight applied to the galaxy. The 
Knax values were ob tained from the NYU-VAGC catalogue (§ 5 of 
iBlanton et alj|2005bh . Figure [2] shows the GSMFs derived from the 



7 Stellar mass estimates from Kauffmann et al. are available from 
http://www.mpa-garching.mpg.de/SDSS/DR4/ 

8 We used the DR4 catalogue of Gallazzi et al. with stellar masses derived 
from z-band Petrosian magnitudes. See footnote |7]f or website. 
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Figure 2. Comparison between galaxy stellar mass functions derived from 
the NYU-VAGC low-2 galaxy sample. The lines represent different meth- 
ods of determining M-L ratios from spectra or colours. The functions were 
determined by summing 1/Vmax in logarithmic mass bins. 



NYU-VAGC using the different stellar mass estimates. For these, 
no weighting was applied to any galaxy (w = 1). 

While all the GSMFs are relatively similar, they all show a 
subtle 'dip' relative to a smooth slope around 10 8 ' 6 .Mo in Fig. [2] 
This feature is a prod uct of using 1 /Knax with large-sc ale struc- 
ture (LSS) variations dEfstathiou. Ellis. & Petersonll988l) . Figure|3] 
shows how the number density of a volume-limited sample of 
galaxies (M r < —19) varies with redshift. There is a clear under- 
dens ity around z ~ 0.013 that is the cause of the GSMF dip. See 
also IBlanton et alj d2005al) . e.g. fig. 9 of that paper, for the effect 
of this LSS on the luminosity function. A simple way of removing 
this effect is to set wi, for each galaxy, equal to 1/71(2) where n is 
the normalised number density as shown in Fig.[3]at the redshift of 
the galaxy. Before recomputing the GSMF in this way, we discuss 
the selection limits related to surface brightness and M-L ratio. The 
mass used for each galaxy in the rest of the paper is obtained from 
the mean of the log M estimates. 



3.1 Selection limits 

One of the key factors in considering the determination of the 
GSMF is the SB limi t, which is often not explicitly identified, 
of a redshift survey dCross & Driven l2002h , From the tests of 
Blant on et all {2005a), as shown in fig. 3 of that paper, the SDSS 
main galaxy sample has 70% or greater completeness in the range 
18-23 mag arcsec -2 for the effective SB fj,RBo,r- 

In order to identify at what point the GSMF becomes incom- 
plete because of the SB limit we computed the bivariate distribu- 
tion in SB versus stellar mass. Figure|4]shows this distribution rep- 
resented by solid and dashed contours (1/Knax and 1/n LSS cor- 
rected). There is a relationship between peak SB and log M, which 
is approximately linear in the range 10 8 - 5 to 10 11 X Q . At lower 
masses, the distribution is clearly affected by the low SB incom- 
pleteness at fj,Rso,r > 23 mag arcsec -2 . Therefore, any GSMF 
values for lower masses should be regarded as lower limits if there 
are no corrections for SB completeness. 

The other important consideration is fact that the r-band selec- 
tion is not identical to the mass selection required for the GSMF. 
This is nominally corrected for by 1/U max but it should be noted 
that galaxies with high M-L ratio at a given mass are viewed over 
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Figure 3. LSS fluctuations in the low-z galaxy sample. The solid line shows 
the relative number density of galaxies brighter than M r = —19 in red- 
shift shells. The number densities were normalised by the volume average 
over z < 0.05. The dashed and dotted lines show the normalised densi- 
ties for M r < —2 0.0 and —21.0, respectively, from the DR4 sample of 
iBaldrvetalJ 12006). The dash-and-dotted lines show the volume averages 
for z < 0.05 (= 1 by design) and z < 0.1 (= 1.13). 
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Figure 4. Bivariate distribution for SB versus mass. The contours represent 
the volume-corrected number densities from the sample: logarithmically 
spaced with four contours per factor of ten. The lowest dashed contour cor- 
responds to 10 — 5 Mpc -3 per 0.5 x 0.5 bin, while the lowest solid contour 
corresponds to 5.6 x 10 -5 Mpc -3 . The grey dashed-line regi ons represent 
areas of low completeness (70% or lower as estimated by Blanton et al. 
l2005ah . The diamonds with error bars represent the median and 1-sigma 
ranges over certain masses with a straight-line fit shown by the middle dot- 
ted line. The outer dotted lines represent ±1 sigma. 



significantly smaller volumes than those with low M-L ratio. Fig- 
ure [5] shows the bivariate distribution in M-L ratio versus mass. 
The limits at various redshifts for the SDSS main galaxy sample 
are also identified. For example, galaxies with M < 10 s Mq and 
\og{M/L r ) > 0.1 are only in the sample at z < 0.008. At these 
low redshifts, the stellar mass and V m &x depend significantly on the 
Hubble-fiow corrections. However, it does appear that the SB limit 
affects the completeness of GSMF values at higher masses than the 
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Figure 5. Bivariate distribution for M-L ratio versus mass. The contours 
represent the volume-corrected number densities: logarithmically spaced 
with four contours per factor of ten. The lowest dashed contour corresponds 
to 10~ 4 Mpc -3 per 0.5 X 0.2 bin. The dotted lines represent the observ- 
able limits for an r < 17.8 magnitude limit and different redshift limits 
(ignoring k-corrections). The grey dashed-line region represents galaxies 
that can only be observed at z < 0.008 where Hubble-flow corrections are 
significant (cz < 2400 km s _1 ). 



M-L limits. At M < 10 8 ' 5 .Mo the SB limit becomes significant, 
while aX M < 10 s Mq, the GSMF is affected both by the con- 
strained volume for high M-L ratio galaxies and more severely by 
the SB incompleteness. 

3.2 Corrected GSMF with lower limits at the faint end 

Figure[6]shows the results of the GSMF determination. The binned 
GSMF is represented by points with Poisson error bars, with lower 
limits represented by arrows. The GSMF has been corrected for 
volume (1/Vmax) and LSS (1/n) effects^ The masses used were 
the average of the four estimated stellar masses. The shaded re- 
gion represents the full range in the GSMF obtained by varying 
the stellar mass used (Fig. [2]l and multiplying the upper number 
densities by 1.13 to account for the expected renormalisation at 
z < 0.1 (Fig.[3j. After renormalisation, the range in f s is 0.038- 
0.048. In addition if we multiply by 1.2 to account for missing light 
(optically-thick regions and corrections to total magnitudes, § [2](, 
we obtain 0.046-0.058. 

The b inned data f or the GSMF can clearly not be represented 
by a single Schechter ( 1976) function. The data were fitted with a 
double Schechter function given by 



• dM 



-M/M* 



AM 

M* 



(2) 



where 4>m AM is the number density of galaxies with mass be- 
tween M and M + dM; with a? < ai so that the second term 
dominates at the faintest magnitudes. Fitting to M > 10 8 .Mq, the 
best-fit parameters are 



9 We compared the GSMF computed using 1/n (2) co rrection for LSS 
variat ions with the stepwise maximum-likelihood method I Efstat hiou et alj 
1988). There was good agreement between the two methods after matching 
normalisations. The former method is simpler when computing bivariate 
distributions. 
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Figure 6. GSMF extending to down to 10 7 Mq determined from the NYU- 
VAGC. The points represent the non-parametric GSMF with Poisson er- 
ror bars; at A4 < 1O 8 5 A^0 the data are shown as lower limits because 
of the SB incompleteness (Fig. |4j- The dashed line represents a double- 
Schechter function extrapolated from a fit to the M > 1O 8 A^0 data 
points. The dotted line shows the same type of function with a faint-end 
slope of a 2 = -1.8 (fitted to M > 10 s5 Mq data). The dash-and-dotted 
line represents a power-law slope of —2.0. The shaded region shows the 
range in the GSMF from varying the stellar mass used and changing the 
redshift range. 



log(M*/M Q ) = 10.648 
(ft/HT 3 Mpc -3 =4.26 
<^/l(T 3 Mpc- 3 = 0.58 



e*i = -0.46 
a 2 = -1.58 



(3) 



with formal errors of 0.013, 0.09, 0.05, 0.07, 0.02. The dashed line 
in Fig.[6]represents this fit extrapolated down to lO 7 ^©. 

Even though the Poisson errors are small, for illustrative pur- 
poses and because systematic errors are clearly significant, we fit- 
ted a function with a 2 = —1.8 fixed. This is represented by the 
dotted line in Fig. [6] and to the eye provides an equally good fit to 
the data at M > W as M & . Given the SB incompleteness (Fig.©, 
a steep faint-end slope such as this cannot be ruled out. 



3.3 Comparison with cluster environments 

The field GSMF shows a clear signal of a change in slope at masses 
lower than the characteristic mass as was alrea dy evident in 
the lum inosity function of the redder SDSS bands telanton et al.l 
l2005ah . Thus there is a significant difference between a faint-end 
slope determined from a Schechter fit around the characteristic 
mass (luminosity) and the faint-end slope at lower masses (lumi- 
nosities). 

Recently, Popes scTet alj d2006t) and ljenkins et al.l ( 120071) have 
confirmed earlier reports of upturns in the faint end of cluster lu- 
minosity functions. These were based on the RASS-SDSS galaxy 



10 It makes minimal difference to the shape of the cosmic volume average 
GSMF, and no difference to the discussion in this paper, if the highest den- 
sity regions (15% of the population, i.e. clusters and compact groups) are 
excluded from the calculation. This justifies the use of the term 'field' to 
describe this GSMF. 
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Figure 7. Comparison of inferred GSMFs from the Co ma Cluster 3.6 -tim 
survey of Ijenkins et all and the SDSS cluster survey of lPopesso et alj with 
the field GSMF (dotted line). The diamonds and trian gles represent the 
Coma core and Coma off-centre columns from table 2 of Jenkins et al, us- 
ing the q — r < 1 .3 restriction for the faint end. The dashed line represents 
the Popesso et al. fit derived from the z-band luminosity function. The data 
were normalised to match the field GSMF around log J\A a = 10.0 ± 0.5. 



cluster survey and 3.6-^im imaging of the Coma Cluster using the 
Infrared Array Camera on the Spitzer Space Telescope, respec- 
tively. 

Stellar M-L ratio variations between cluster galaxies are typ- 
ically less severe than between field galaxies. Using a simple con- 
version between absolute magnitude and stellar mass given by 
log Ms = (M so i ar - M)/2.5 + \og(M s /L) , we converted the 
cluster luminosity functions to GSMFs: with M ZjSO i ar = 4.4 (AB 
mags), \og(M a /L z ) = 0.2 (solar units), M3. 6iSO i a r = 3.3 (Vega 
mags) and \og{M a / L^s) = —0.5 (solar units). The conversion 
factors were estimated using PEGASE and the filter curves. Figure|7] 
shows the resulting cluster GSMFs. The faint-end upturn is evident 
at .M s < 1O 9 X and is significantly steeper than the field GSMF. 
Thus it appears that the slope of the GSMF around 10 8 to IQ^Mq 
depends significantly on the environment. Note however that these 
cluster results rely on estimated subtraction of background galaxy 
counts. 

As well as the difference in the slope of the GSMF in the range 
10 s to 1O 9 X , there is the more established difference between 
the morphologies of these low-mass galaxies in different environ- 
ments. In clusters, these are predominantly dwarf ellipticals (dE); 
whereas in the field, these are predominantly late-type spirals (Sd) 
and irregulars (Sm,Im) (e.g. by estimatin g stellar masses for galax- 
ies in the Nearby Field Galaxy Survey of ljansen et alj |2000). Field 
low-mass galaxie s are generally form ing stars and have substantial 
reservoirs of gas ( Sw aters et al.ll2002h . and therefore their baryonic 
masses can be several times their stellar masses. 



4 THE STELLAR MASS-METALLICITY RELATION 
AND THE BARYONIC MASS FUNCTION 

In order to convert our stellar mass function (MF) to the more fun- 
damental baryonic MF, we develop a method for deriving the stellar 
mass fraction (i.e., conversion factor of baryonic mass to stars) in 
terms of the stellar mass. This can be achieved by using the well 
established relation between stellar mass and metallicity coupled 
with a metallicity to stellar mass fraction relation, which can be de- 
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termined from a simple model. This method is laid out below. We 
assume the following: (i) gas in a galaxy is well mixed, in partic- 
ular, we do not include metal-enriched outflows in the model; (ii) 
the galaxy gas-phase metallicity measured from an SDSS spectrum 
represents an effective average over the whole galaxy and can be 
related to the global gas fraction, i.e., there is no consideration of 
metallicity and related gradients in galaxies. Despite these simpli- 
fications, the derived average stellar mass fractions are shown to be 
consistent with direct measures of gas masses. 

4.1 Relating metallicity to stellar mass fraction 

The importance of stellar mass is demonstrated by the tight 
relationship between gas-phase metallicity and stellar mass 
dTremonti et al.l2004h . and the relationship between metallicity and 
the fraction of baryonic mass in stars. The latter is illustrated us- 
ing a three-component model with stellar mass, retained gas and 
expelled gas, which reduces to the closed-box model in the case 
of zero expelled gas. Regardless of the problems with this simple 
model, there clearly must be a fundamental relationship between 
the fraction of baryons locked up in stars, gas and the progress of 
chemical evolution in a galaxy . 

From eq. 5 of Edmunds (1990) with a simple outflow that is 
proportional to the SF rate and no infall, we can set 

d(gZ) = [ap -aZ - oZ) ds (4) 
dg = — (a + o)ds (5) 

where: g is the mass of gas, Z is the metallicity of the gas, p is 
the yield defined as the mass of metals produced per mass of long- 
lived stars formed, s is the integrated mass of stars formed, a is the 
fraction of mass in stars that is not instantly recycled, i.e. as is the 
stellar mass of the galaxy, and o is the mass of gas expelled per 
mass of stars formed. The integral is then given by 



1.00 



dZ 
Z 



\a + oJ 

Win 

\a + oJ 



dg 



(6) 
(7) 



where gi is the initial gas mass, i.e., the total mass of stars and gas 
(the baryonic mass), and g = g\ — as — os . Setting o = 0, Eq.[7j 
reduces to the standard closed-box solution given by 



pin - 

g 



-pin (1-22 



(8) 



Figure [8] shows how the fraction of mass in stars (as/gi), re- 
tained inter-stellar gas (g/gt) and expelled gas (os/gi) depends on 
the metallicity. The curves were obtained from Eq.UJwith a — 0.6 
and o = 0.0, 0.4, 1.2, 3.0. The curves are shown for g/gi ^ 0.01 
i.e. a minimum of 1% of the mass remaining as inter-stellar gas. 

The simple model demonstrates that the gas-phase metallicity 
can be a relatively accurate predictor of a galaxy's SF efficiency 
(as/gi or baryonic-to-stellar mass conversion factor) regardless of 
a wide range of outflow scenarios. If the retained gas is greater in 
mass than the expelled gas, the SF efficiency will be within 30% of 
the closed-box estimate derived from the metallicity regardless of 
the outflow factor. If there is significant expelled gas, using Eq.[8]to 
relate metallicity to the stellar mass fraction will result in an under- 
estimate of as /(as + g) and an overe stimate of as /eg . If there ar e 
continuing infalls and outflows (e.g. |Palcanton]|2007l ; lErbll2008l) , 
the case is less clear because of the dependence on the metallicity 
of the infalling gas. 
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Figure 8. The mass fraction of stars (solid lines), retained inter-stellar gas 
(dashed lines) and expelled gas (dotted lines) versus the gas-phase metal- 
licity for the three-component model (Eq.|7J. The darkest lines represent a 
closed-box model with no outflow while the lighter lines represent increas- 
ing outflow factors (0, 0.4, 1.2, 3.0). The diamonds and squares represent 
points at which the mass in expelled gas is equal to the mass in retained gas. 
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Figure 9. The s tellar mass-metallicity relation derived from the emission- 
line analysis by[jremonti et al. 12004). The squares represent the median 
while the lines represent the 16th and 84th percentiles in logarithmic mass 
bins. The dotted lines sh ow the solar abundance value determined by 
lAllende Prieto et all i200lll . 



Figure [9] shows stellar mass-metallicity relation using 
gas-ph ase oxygen abundances estimated by iTremonti et al.1 
( 2004) . This does not vary signi ficantly with environment 
dMouhcine, Baldrv. & Bamfordl 2007) suggesting that there is a 
fundamental relationship bet ween the present-day stellar mass of 
a galaxy and its SF efficiency. iBrooks et al.1 12007) concluded from 
their modelling that "low star formation efficiencies ... are primar- 
ily responsible for the lower metallicities of low-mass galaxies and 
the overall M-Z trend". 

Note that the estimated metallicit y can depend on the emission 
lines considered and the calibration |K obulnickv & K ewley||2004l : 
ISavaglio et alj|2005l ; iKewlev & Ellisorj|2008h . However our con- 
sideration here is only that the measured M-Z relation implies an 
increase of average gas-phase metallicity with mass, and that the 
relation between metallicity and stellar mass is fairly tight. 



© 2008 RAS, MNRAS OOO.miBl 



8 /. K. Baldry et al. 



4.2 Relating stellar mass to baryonic mass 

The variation in SF efficiency implied by the M-Z relation allows 
one to estimate a baryonic MR We assumed (i) the median mea- 
sured oxygen abundances at a given stellar mass are representative 
of the average gas-phase metallicity in the inter-stellar medium, (ii) 
the yield is independent of galaxy mass and time, (iii) the closed- 
box model can be used to relate SF efficiency to oxygen abundance, 
and (iv) the SF efficiencies of the most massive galaxies are about 
90%. Next, we defined a parametrisation to relate the total baryonic 
mass of a galaxy (Mb) to the SF efficiency: 



Ms 
M b 



M b 



Mb + Mo 



(e„ - ei) + e t 



(9) 



where Mo is a characteristic mass, and e u and ei are upper and 
lower limits of the SF efficiency. Combining with Eq. [8] this 
parametrisation was fitted to the M 3 -Z data (using O/H, Fig. [9] for 
Z). Figure [Tol shows the SF efficiency versus log Ms and log M b 
from the best-fit parametrisation (log.M = 9.6, e u = 0.90, 
e; = 0.09, 12 + logp = 8.84 where p is the oxygen yield). 

The implied relationship between stellar and baryonic mass 
was then used to compute the baryonic MF (from the GSMF or 
from the stellar masses of all the galaxies). This is shown by the 
dashed line in the lower panel of Fig. \10\ The solid line is a fit to 
the MF at Mb > 10 8 5 Mq. The faint-end slope in this case is 
given by a steep Q2 — —1.9 (parameters as per Eq.|3]are 10.675, 
4.90, -0.42, 0.61, -1.87.) 

The steep part of the baryonic MF is caused by the initially 
fairly steep GSMF at Ms < W 9 Mq coupled with an increase in 
gas-to-stellar (G-S) mass ratio from ~ 1.5 at 10 9 .Mq to ~ 4 at 
10 8 X Q to ~ 8 at 10 7 Mq (Fig.QJJ dotted line in upper panel). 
Though the G-S mass ratios for low-mass galaxies are large they 
are in fact similar to those obtained directly from HI surveys. For 
example, the abov e values are similar to the ratios estimated by 
iKannappanl d2004l) for the low-mass end of the blue sequence. 

In order to test the average G-S mass ratios implied by 
our parametrisation, we used stellar a nd atomic gas masses de- 
rived from the Wester bork HI Survey dSwaters & Balcells |2002|; 
iNoordermeer et al"]|2005l) and the literature compilation of Garnettl 
(2002). The stellar masses were estimated using the simple rela- 
tio ns of | Bell et al] (1200 3b) applied to the B and R ph otometry 
of ISwaters & Balcells! and the B and V photometry of iGarneut 
and using \og(M s / 'Lb) = 0.3 for the early-type spirals of 
INoordermeer et aU We al so matched the HI Parkes All-Sky S ur- 
vey (HIPASS) catalogue dMever et alj|2004 IWong et alj|2006h to 
the NYU-VAGC catalogue. This resulted in 170 galaxies with one- 
and-only-one match within 0.25° and Azc < 250 tons" . 

Figure [TTJ shows average G-S mass ratios in bins of stellar 
mass for these surveys; also shown is the relation derived from our 
parametrisation. The HIPASS mass ratios lie above this line, not 
surprisingly, because the HI selection misses galaxies with low G- 
S mass ratios. Thus these derived mass ratios should be regarded 
as upper limits. The G-S mass ratios derived from the optically- 
selected Westerbork survey are in good agreement, which lends 
support to our parametrisation of the SF efficiency. The Garnett 
compilation points to a flatter relation but this was using a signifi- 
cantly smaller, heterogeneously-selected sample. 

Various authors have found that the effective yield is 
lower in low mass galaxies dGarnettll2002l ; IPilvugin et alj|2004 
iTremonti et alj|2004h with an interpretation bein g that there is sis 
nificant metal loss by metal-enriched outflows i Dalcantorj 12007 
The effective yield is defined as 
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Figure 10. SF efficiency versus mass and the implied baryonic mass func- 
tion. Upper panel: The model SF efficiency is defined as the mass in stars 
formed divided by the mass available to form stars over the life of a galaxy. 
This is plotted versus stellar mass (dotted line) and baryonic mass (solid 
line). The squares represent the O/H values of Fig. [9] converted using a 
monotonic relationship to efficiency determined from a closed-box model 
and plotted against stellar mass. Lower panel: The dotted line represents the 
GSMF from Fig. [6] The dashed line represents the implied baryonic mass 
function given the relationship between SF efficiency and stellar mass. The 
solid line represents a double Schechter fit to these data at M b > 10 9 Mq 
with faint-end slope ct2 — —1.9. The dash-and-dotted line is this fit con- 
verted back to stellar mass. The arrows illustrate the relationship between 
baryonic and stellar mass in this simple model. 



p e ff = Zj In 



as + g 



(10) 



such that p c ff = p for the closed-box model, and p c g < p with out- 
flows. Our analysis uses SDSS-aperture metallicities and a closed- 
box model to predict average gas masses. This gives approximately 
correct global G-S mass ratios at least in comparison with the 
Westerbork samples (Fig. [Til , which is consistent with the vari- 
ation in SF efficiency wi th mass being the primary cause of the 
Ms-Z relation (see also iBrooks et al J 120071 : lEitison etai]|2008l: 
iMouhcine et"ai]|2008l : Pfassis et alj|2008l) . This does not mean that 
there are no significant outflows, only that they may be a secondary 
cause of the Ms-Z relation. 

The precise shape and steepness of the baryonic MF does de- 
pend on the detail of the G-S mass ratios and the dispersion in the 
relationship with stellar mass. However, our main aim is to illus- 
trate the relationship between galaxy baryonic mass function and 
GSMF, and so we use the simple stellar-to-baryonic mass relation 
derived here. We also note that while the change in slope convert- 
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Figure 11. Average atomic G-S mass ratios. The symbols represent G-S 
mass ratios from Hi measurements with a 1.33 correction factor to account 
for Helium. The data from 4—25 galaxies were averaged for each point. The 
vertical error bars represent the standard errors; while the horizontal error 
bars represent the ranges of each sample. See text for further details. The 
dashed line represents the parametrisation of Eq.[9]assuming the closed-box 
model and that all the gas is atomic gas. The dash-and-dotted line represents 
the same parametrisation but assuming 70% of the gas is atomic gas. Note 
that gas-poor ellipticals are not considered in these samples, however, they 
are only a significant field population at > 10 11 Mq and so their gas prop- 
erties can be ignored in terms of computing the baryonic MR 



ing from the GSMF to the baryonic MF could be exaggerated, the 
GSMF slope may be underestimated and so the faint-end slope of 
the baryonic MF could be ~ —1.9 even accounting for a more 
complicated relation. 



4.3 Comparison between galaxy and halo baryonic mass 
functions 

The implied baryonic mass function (Fig. 1 10b can be considered 
to be the sum of stars and any gas involved with the cycle of SF, 
in and around each galaxy, as long as the relationship between 
Ms I Mb and Z is not strongly affected by o utflows. Other esti - 
mates of the baryonic M F ha ve been made by Bell et al] d2003al) . 
iRead & Trenthaml d2005h and lShankar et al.l ( 2006). These include 
stars and atomic gas in galaxies, and molecular gas in the case of 
the first two estimates. 

Figure [T2] shows a comparison between th e implied baryoni c 
mass function and the diet Salpeter version from lBell et afld2003al) . 
There is a significant offset between them, which can be reconciled 
by (i) allowing for the number density and missing light corrections 
of § !3.2l and (ii) slightly lowering the masses of lBell et alj because 
the stellar M-L ratios appear on the high side (§|2}. This is partly be- 
cause of using the diet Salpeter IMF compared to those of Kroupa 
or Chabrier (Appendix). After plausible corrections, the galaxy 
baryonic MFs are in good agreement. There is no clear evidence 
for a steep slope at M < 10 9 ' 5 .M© in the iBell etaU function, 



but given the lar ge error bars, it is consis ten t with the steep slope 
we find. Neither lRead & Trentharr] d2005t) or lShankar et alj'd2006h 
noted such a steep slope in the baryonic MF. IShankar et all 's G-S 
mass ratios were based on a calibrat ion of HI and stellar ma sses as 
a function B-band luminosity (from lSalucci & Pe rsic 1999). Their 
average G-S mass ratio is ~ 1 at Ms = 10 8 .Mo, which is signifi- 
cantly lower than our estimate. 

For comparison, halo MFs are shown in Fig. [T2] assuming 



Figure 12. Comparison between galaxy and halo baryo nic mass functions . 
The thick dotted line represents the MF determined bv lBell et alj d2003at) ; 
with the thin dotted line representing masses reduced by 15%. Error bars 
are shown for the low-mass end. The thick solid line represents the MF 
from Fig. 1101 with the thin line representing masses increased by 20% 
and number densities by 13% (cf. § 13.2) . Halo MFs with masses multi- 
plied by fli,/Q m are shown by the dashed lines (Millennium Simulation; 
IShefh&Tormenll999l using longer dashes), and by the dash-and-dotted line 
for the 'galactic halo MF' IShankar et al.l2006h . 



a constant fraction of mass in baryons in each halo (f2(,/O m ). 
The long-dashed line was derived from the simulations of 
ISheth&Torrnenldl999l) while the dashed line was derived from the 
Millennium Simulation dSpringel et al.ll2005l ) F^l The SF efficiency 
implied by the M a -Z relation shows that the galaxy baryonic MF 
is approximately as steep, at Mb < 10 9 ' 5 Mp, as the halo M F. 
The 'galactic halo MF' given by eq. 9 of IShankar et alj d2006t) is 
also shown. This galactic halo MF, which includes halos and sub- 
halos hosting a galaxy but not group and cluster halos, is very sim- 
ilar in shape to the standard halo MF at Mb < W^Mq. This is 
because the slopes of sub-halo MFs do not depend strongly on the 
mas s of the main halo and have values sim ilar to that of the halo 
MF dPe Lucia et al.ll2004l;lRe~ed et al.ll2005l) . 

To illustrate the implications of the form of the mass func- 
tions on galaxy formation efficiency, we computed the efficiency 
as a function of halo baryonic mass required to reproduce the 
galaxy MFs. This assumes a one-to-one and monotonic relation- 
shi p between halo mass and galaxy mass. In detail this is described 
by IShankar et al. 1 20061). a nd ea rlier using galaxy lumi nosity by 
iMarinoni & Hudsonl ( I2062I) and IVale & Ostrikeil d2004l) . For the 
halo MF, we use the galactic halo MF with mass multiplied by 
fii,/fi m . Figure [T3l shows the efficiency versus halo mass and the 
reconstructed galaxy MFs. The efficiency represented by the solid 
line can be regarded as the fraction of baryons that fall onto a galaxy 
and are available for forming stars; while the dashed line represents 
the fraction of baryons that have formed stars. These reproduce the 
galaxy baryonic MF and GSMF, respectively. 

The figure demonstrates the implication that there is a level- 
ling off of the baryonic-infall efficiency at low masses, while SF 
efficiency continues to fall, and shows a peak galax y formation ef- 
ficienc y at 10 U1 M S (1O 119 A4 including CDM). IShankar et al.1 
(2006) also looked at SF efficiency versus halo mass and demon- 



The Millennium Simulation is available at 
http://www.mpa-garching.mpg.de/millennium/ 
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Figure 13. Idealised reconstruction of galaxy formation efficiency. The top 
panel shows the fraction of baryons falling onto a galaxy (solid line) and the 
fraction of baryons forming stars (dashed line). The lower panel shows the 
result of using the efficiency curves to go from the galactic halo baryonic 
MF (dash-and-dotted line) to the galaxy baryonic MF (solid line) and to the 
GSMF (dashed line). A standard halo baryonic MF that does not distinguish 
between galactic, group and cluster halos is shown by the dotted line. The 
arrows illustrate the relationships between the mass functions. They are not 
flat because the number density per logarithmic bin changes depending on 
the derivative of efficiency with respect to mass. 

strated a similar result (their fig. 5) albeit with a flatter falloff at 
higher masses and a steeper falloff at lower masses related to the 
different GSMF utilised by them. The reduced efficiency above 
and below this scale are inherent in recent semi-analytical models 
of gal axy formation including SF and active galactic nuclei feed - 
back teower et al. 2006; Cattaneo et al. 200rj: lcroton etai] |2006). 
The scale is related to a minimum in feedback effects (cf. fig. 8 of 
iDekel & Birnboimll2006l) . 



5 SUMMARY AND CONCLUSIONS 

The field low-redshift GSMF has been generally measured down to 
~ 10 8 ' 5 Mq (Fig. [I). This accounts for the majority of stellar mass 
in the universe. The SMD of the universe is in the range 5-7 per 
cent of the baryon density (conservatively 4-8 per cent, assuming 
SI;, = 0.045). The sources of uncertainty are corrections for dust 
attenuation, underlying stellar M-L ratios, missing low-SB light, 
and normalisation of the GSMF. 

We determined the field GSMF using the SDSS NYU-VAGC 
low-redshift galaxy sample. At masses below ~ 10 9 Mq, the shape 
of the GSMF computed using a standard 1 / V max approach was af- 
fected by LSS variations. This was corrected for using number den- 
sity as a function of redshift of a volume-limited sample of galaxies 
(Fig. ©. Analysis of the bivariate distribution of stellar mass ver- 



sus SB indicated that the number densities should be regarded as 
lower limits at masses below ~ 10 8 ' 5 Mq (Fig. 0. Despite this 
incompleteness, there is clear evidence for an upturn in the GSMF 
(Fig. [6} with a faint-end slope 02 — —1.6 (Eq.[3](. This represents 
the power-law slope over the range < 10 8 to ~ 10 9 .M©. Steeper 
slopes are also plausible. Slopes of ~ — 2 have been measured in 
clusters over the same mass range (Fig. [7}. The processes shap- 
ing the GSMF in clusters will include stripping of gas and stars, 
while feedback processes may be more dominant in shaping the 
field GSMF. 

At masses below 10 10 Mq, blue-sequence (late-type) galaxie s 
are the dominant field population (e.g. fig. 10 of lBaldrv et al.l2004h . 
The gas-phase metallicity of these galaxies can be related to a SF 
efficiency that is defined as the total mass of stars formed divided 
by the available baryonic mass for forming stars. The closed-box 
formula can remain a good estimate of the SF efficiency even with 
moderate outflows (Fig. |8j, e xcluding the cas e of metal-enriched 
outflows for gas-rich systems (Dalcanton 2007). 

Using a simple relationship between stellar mass and bary- 
onic mass, based on the M s -Z relation and the closed-box formula 
for Z that neglects outflows, we converted the field GSMF to an 
implied baryonic mass function (Fig. |10t . The resulting faint-end 
slope 0:2 ~ — 1.9 is similar to the halo MF. We note that this is only 
suggestive as it depends on the form of the conversion between M s 
and Mb and the dispersion in this relationship (Fig. II It. The shape 
of the galaxy ba r yonic M F is consistent with the non-parametric 
MF of lBell et alj i2003al) (Fig.[l2j. 

Taking the shape of the implied baryonic mass function at face 
value, we compared this with a simulated halo baryonic MF. Using 
a one-to-one relationship between halos and galaxies, these can be 
used to determine the galaxy formation efficiency (the fraction of 
baryons falling onto a galaxy) as a function of halo mass (Fig.|13ll. 
This illustrates how varying efficiency with mass can be used to 
obtain galaxy mass functions from the halo MF, or from a similarly- 
shaped sub-halo MF. The peak in the formation efficiency curve 
may correspond to a minimum in feedback efficiency. 

There is no evidence yet of any cutoff in mass, below whic h 
baryons do not collapse into galaxies (cf. IDekel & Wool 120031) . 
Rather we find the baryonic-infall efficiency levels off to ~ 10% 
rather than continuing to plummet with mass. It is possible that a 
cutoff mass scale, imprinted in the shape of the field galaxy bary- 
onic MF, could be found by future deeper surveys. To robustly iden- 
tify this scale requires: (i) wide-field deep imaging with reliable 
identification of galaxies down to at least fi r ~ 25 mag arcsec 2 ; 
(ii) spectroscopy of large samples down to r ~ 20 mag, or a simi- 
lar effective depth in near-IR selection, in order to obtain accurate 
distances (z > 0.01) and metallicities for low-mass galaxies; (iii) 
wide-field HI surveys in order to estimate gas masses more directly. 
The prospect of such a measurement within the next decade is good 
with the advent of new wide-field optical instruments and survey- 
efficient radio telescopes. 
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APPENDIX A: MASS-TO-LIGHT RATIOS OF STELLAR 
POPULATIONS 

Stellar M/L ratios of galaxies are generally based on evolution- 
ary population synthesis with the basic ingredient being sim- 
ple stellar populations (SSPs). This section outlines some is- 
sues related to M/L determination: defini tion, IMF, metallicity , 
age, bursts, dust, synthesis c ode. [See also IPortinari et all (2004): 
ItCannappan & Gawised d2007h .1 

The following are common assumptions for estimating stellar 
masses. The IMF is valid from 0.1 to 100 (or 12 0) Substel - 
lar objects, < O.IMq, are not included; from the lChabrierl (2003) 
IMF, these add up to 5-10% of the stellar mass. Stellar remnants 
(white dwarfs, neutron stars, black holes) are included; this is typ- 
ically a 10-20% factor in the PEGASE models. Stellar mass is the 
remaining mass in stars and remnants as opposed to the integral of 
the SF rate, i.e, M. = (1 — R) J SFR where R is the recycled frac- 
tion. This is the definition used here and is appropriate considering 
the analysis of §|4] 

A significant factor is the choice of IMF. Figure lATl shows how 
M-L ratios in the r- and J\"-bands vary with IMF choice. The diet 
Salpeter has the highest masses considered here as it was calibrated 
to have ma ximum M-L ratios c onsistent with dynamical mass 
constraints teell&de Jon3l200ll) . Thus on the assumption of a 
universally-applicable IMF, this is the most-massive IMF at a given 
SSP colour that is valid. At the high-mass end of the IMF, steep 
slopes implying fewer high-mass sta rs are ruled out by measure - 
ments of cosmic luminosity densities ( iBaldrv & Glazebrookl2003h . 
Thus reasonable M-L ratios are obtained with IMFs of Kroupa, 
Chabrier (0.91-0.93), diet Salpeter (1.07-1.11), Baldry & Glaze- 
brook (0.90-1.04), and Kennicutt (0.74-0.88), where the ranges in 
parentheses are derived stellar masses relative to the Kroupa IMF 
at SSP colours evaluated over the range of 0.2-0.85 in q — r. [If 
the IM F varied w ith galaxy luminosity (Hoversten & Glazebrookl 
l2008h . over time jDave1l2008l ; Ivan Dokkumll2008l) and/or in star 
bursts dFardal et al.ll2007l) . this would clearly complicate the deter- 
mination of stellar masses.] 

A more significant consideration is the choice of prior (al- 
lowed SF histories, etc.) and population synthesis code. The follow- 
ing figures are provided to illuminate some of these considerations. 
Figure lA2l shows M-L ratio versus co lour for SSPs over a range 
of metallicit ies de rived from PEGASE dFioc & Rocca-Volmerangd 
1 19971 1 19991) and iBruzual & Charlotl d2003h models. Figure |A3] 
shows the effect of adding a 100-Myr burst contributing 5% of 



the stellar mass. Figure lA4l shows a comparison between PEGASE 
and preliminary M-L ratios derived from BaSTI (M. Salaris, private 
communication; IPietrinferni et al.l2004l) . Figure |A5l shows PEGASE 
models for SF histories with constant rate. 

Last but not least, there is the complication of dust attenuation. 
The arrows in Figs. IA2HA5I show the effect of 1 magnitude of at- 
tenuation, at the V band, on the colours and M-L ratio s: both for a 
Small Mag ellanic Cloud (SMC) screen law dPeil 1992b and a A" 7 
power law dCharlot & Fall 2000). The latter law allows for a range 
of attenuation to different parts of a galaxy making it 'greyer' than 
a screen law. 

An example of the importance of the prior and dust law is 
demonstrated by the fitting to the photometry of the NYU-VAGC 
sample (§ [3}- Fitting only to photometry, as opposed to spectral 
features, has the advantage of being less sensitive to aperture bias 
but the disadvantage that colours are highly sensitive to dust atten- 
uation. Various fits including different dust laws were tested (also 
fitting to ugriz + JHK). It was found that using a screen law sig- 
nificantly lowered the mass of luminous red galaxies (early types), 
which were fitted with significant dust attenuation, compared to 
assuming no dust. The lower masses occurred because younger- 
age models with dust reddening to reproduce the observed colours 
give lower M-L ratios than older-age models with no dust, i.e., the 
screen-law dust vector (change in M-L ratio as a function of colour) 
is shallower than the age vector for red galaxies (Fig. IA2t . Using 
a A -0,7 dust law did not lower the masses but early- type galax- 
ies were still fitted with significant dust, which is not consistent 
with our general knowledge of these galaxies, unless the fitting was 
restricted to solar metallicity only. Comparable GSMFs were ob- 
tained either by allowing for varying dust attenuation and assuming 
solar metallicity or by allowing for a range of metallicities and as- 
suming no dust. The former is more appropriate for the high-mass 
galaxies while the latter is more appropriate for low-mass, and low- 
metallicity, galaxies. 
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Figure Al. M-L ratios of simple stellar populations computed using PEGASE with ages from 100 Myr to 12Gyr versus colour at solar metallicity. The 
left panel shows r-b a nd M-L ratios while the right panel shows K -bmvd M-L ratios. The tracks represent different IMFs jKroupjliOQll; IChabrieJl200i 
iBell & de Jon3l200ll : liaidrv & Glaze brook 2003; Kennicutl| [l983h . The upper-mass slopes of the IMFs are T = 1.3, 1.3, 1.35, 1.15, 1.5, respectively 
(dn/dlog m oc m~ r ). The low-mass slopes are reduced at m < 0.5 or < 1A4 q except for the diet Salpeter, which uses a 0.7 correction factor. 




Figure A2. M-L ratios of simple stellar populations with ages from 100 Myr to 12Gyr versus colour. The left panel shows r-band M-L ratios while the right 
panel shows i\~-band M-L ratios. The tracks represent different metallicities colour-coded according to the legend in the top-left of each panel. The arrows 
correspond to the effect of A v = 1 mag of dust attenuation for an SMC screen law (thicker line) and a A - 0,7 power law. 
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Figure A3. As per Fig. lA2l except a 'burst' corresponding to an SSP of 100 Myr contributing 5% of the stellar mass has been added to each population. Since 
the SSPs are plotted from 100 Myr, the tracks start in the same position as Fig. lA2l 
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Figure A4. As per Fig. lA2l except comparing M-L ratios from BaSTI with PEGASE. The BaSTI models cover a larger metallicity range than implied by the 
legend with [Fe/H] from 0.4 to -2.3. 
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Figure A5. As per Fig. lA2l except for constant-rate SF histories. 
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